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Abstract 

The gravity water wave black (GWBH) hole analog discovered by 
Schutzhold and Unruh (SU) is extended to allow for the presence of 
turbulent shear flow. The Riemannian geometry of turbulent black 
holes (BH) analogs in water waves is computed in the case of lami- 
nar tirbulent shear flow. The Riemann curvature is constant and the 
geodesic deviation equation shows that the curvature acts locally as 
a diverging lens and the stream lines on opposite sides of the analog 
black hole flow apart from each other. In this case it is shown that 
the curvature quantities can be expressed in terms of the Newtonian 
gravitational constant in the ergoregion. The dispersion relation is 
obtained for the case of constant flow injection. 



1 Departamento de Ffsica Teorica - Institute- de Fi'sica Rua Sao Fco. Xavier 524, Rio 
de Janeiro, RJ 

Maracana, CEP:20550-003 , Brasil. 
E-Mail.: garcia@dft.if.uerj.br 



1 Introduction 



Schutzhold and Unruh [1] have built water gravity wave analogs of black 
holes. Since the gravity water waves speed can be adjusted easily by varying 
the height of the fluid, these scenarios possess certain advantages over the 
more well known dielectric and acoustic black holes analogs. In acoustic 
BHs the acoustic horizon appears if the velocity of the fluid exceeds the 
speed of sound in the medium of the liquid. The Riemannian geometry of 
vortex acoustics in such vortex sinks [2] has been investigated in detailed by 
Fischer and Visser [3, 4]. Recently extension of this geometry to the non- 
Riemannian geometry of vortex acoustics to taking into acount rotating non- 
viscous and rotating viscous fluids have been discussed by Garcia de Andrade 
[5, 6]. In this paper we investigate the detailed geometry of the turbulent 
GWBH effective metric. As an example we apply Riemannian geometry to 
the turbulent acoustic shear flow recently investigated by J. de Barros in her 
DSc. thesis [7]. Her model is composed of shear laminar flows in between two 
other fluid layers. The shear layers here represents the surface of the water 
from above and the bottom of the water tank on the other side of the fluid. 
The paper is organised as follows: In the section two we present the physics 
of the model [7] as well the metric of the turbulent fluid shear layer endowed 
with an ergoregion of the analog GWBH. In section 3 the computation of 
Christoffel connections and Riemann and Einstein tensors are given and the 
application to a shear laminar turbulent flow is given which allows us to 
compute the curvature quantities in terms of the Newtonian gravitational 
constant. In section four the computation of the dispersive relation of the 
turbulent irrotationally perturbed flow in the case of stationary injection 
of mass into the fluid, and the group velocity of the flow. Finally in the 
conclusions we discuss applications and future prospects. 

2 The turbulent fluid model and the water 
GWBH 

Turbulent phenomena are in general quite complicated and can only be 
treated numerically and seldom treated analytically [8]. Nevertheless some 
interesting particular models of shear flow turbulence can be treated analyt- 
ically such as the case of laminar shear flows in between two layers where 
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the background turbulence is perturbed. In this section we shall consider 
one of these models described by J.D. Barros [8] and show that they can be 
given a geometrical interpretation as the Riemannian geometry of a GWBH. 
The dynamical equations are given by the modified Euler and conservation 
equations 

^ + pV.v = q(t) (1) 
P^ t + VP = / (2) 

— * 

where the source / is a dipolar source while the turbulence injection is given 
by the function q(t). The differential operator is given by 

§ = (| + ff(, llr ).V) (3) 
Linearization of the dynamical equations of the flow according to the scheme 

P = Po + P (4) 

p = Po+p' (5) 

v = U + v> (6) 

where Vpo = = Vpo, where cq is the speed of sound in the medium and 
U = U(x,y)e[ is the velocity of the background turbulent flow, yields the 
following dynamical flow equations 

% + A)V.t? = q (7) 

Po(^) + (VU.v')ei + W = f (8) 

To simplify matters we shall assume that the flow is irrotational or v' = Vip 
where ip is a scalar function. This assumption is reasonable since we shall 
next consider a shear thin layer where no bulk vorticity exists and only 
localized vortices are appreciable. We also assume that the external force 
can be obtained by the relation / = V0 and VU.v' = 0. This term can 
be incorporated again when we investigate the non-Riemannian geometry 
of turbulent flows [3] later on. The external force here is similar to what 
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happens in the Newtonian gravitational analogue of gravity waves [4] . After 
some usual algebraic manipulations one obtains 



By making use of expression (3) into (1) yields the generalised wave equation 
for the turbulent flow 



c 2 V 2 ^ -d + U.V)^ = C -[q- ^ A] (10) 



'dt ' " " ' T po 1 * 1 Dt K c 2 

Throughout the calculations we also use the following relation 

D lP ' _ 1 D lP ' 
~Dt ~ ^~Dt 



(11) 



This covariant scalar wave equation (10) gives rise to the following effective 
metric 

,00 



9 U \ff = 1 (12) 

9°\ff = V (13) 

and 

g v eff = U l W - c 2 5 lj (14) 

where 5^ is the Kronecker delta and = 1,2,3). The line element repre- 
senting the analog effective metric is 

U 2 

ds 2 = -(1 - —)dt 2 - 2Uidx l dt + Sijdx'dx 3 (15) 

This "spacetime" region represents the PGL type metric of the analog BH 
which reduces to the Schutzhold-Unruh metric of gravity water waves when 
c 2 = gh where g is the usual Newtonian gravity acceleration and h is the 
height of the tank of water. From the above metric one notes that the 
ergoregion of the effective BH is given by U 2 = c 2 . 
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3 Riemannian curvature of the GWBH and 
gravitational analog of diverging lens 

In this section we shall compute following reference [9] the curvature quan- 
tities associated with the GWBH. We also compute the geodesic deviation 
equation and show that the GWBH curvature acts locally as a diverging lens 
with a complex curvature frequency. Let us now consider the Riemann cur- 
vature components by applying the Gauss-Codazzi decomposition used by 
Fisher and Visser in reference [9]. Applying their expressions to our partic- 
ular example one obtains, first the Riemann- Christoffel connections 

Y t ij = D i:j = (diUj + djUi) (16) 

Y\ t = U t U 3 D l3 = \{U.V)U 2 (17) 

T'u = -UjDij (18) 
T\ t = -UjDij (19) 
r jk = U t D jk (20) 

T\ t = -m - \{5 l3 - ITU^djtJ 2 (21) 

where the first term vanishes since U l does not depend on time 

V tj = -U l U k D jk (22) 

where we have put the vorticity £1^ to vanish since there is no vorticity in the 

— * 

Lily model considered here and the perturbation v' is irrotational. With these 
tools in hand we are able to compute the Riemann curvature components for 
the model where we consider the U = U(x,y)ei where U = By , B being a 
proportionality constant, one obtains the following curvature components 

Rijki = {DkiDij — DjiDik) (23) 

Rtijk = UiiDkiDij - D 3l D ik ) (24) 
Rtuj = -{D% - U k Uk,ij + UMDkiDij - DjkDu) (25) 

where we have used also the expression TrD = divU = B. To simplify the 
computations we consider the Riemann components in a orthonormal tetrad 
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Minkowski frame (e"), where fi — 0, 1, 2, 3. The metric in the tetrad frame is 
rjab where (a = 0, 1, 2, 3), is given in terms of coordinate Riemann metric as 

9^v = 7?afce%e^ ( 26 ) 

where we use the same simple gauge used by Fischer and Visser [9] as e*' t = 
— 1, e l t = —v l and e J j = 8^ Finally the Riemann components in this gauge 
are 

Rfy't'y' = —D yx D xy = —B 2 (27) 

Rx'y'x'y' = —D 2 yx = ~R 2 (28) 

These spacetime reveal that this particular gravity wave black hole analog 
with laminar shear flow has a constant Riemann curvature. In this same 
gauge the Ricci components are 

Re* = -V 2 U t = -^(Vx^ = (29) 
and the only non-vanishing components of the Einstein tensor are 

GW = ~B 2 (30) 

G x >y> = B 2 (31) 

With these curvature components we are now able to calculate the geodesic 
deviation equation 

n 2 

U - - - t'\2i 



d\ 2 



[5y] + [(Ry't'y't' + Ry'x'y'x'W Y\[8y\ = (32) 



where 5y is a space-like separation of a family of geodesies in the x — direction 
with tangent vector u a = («*',«*', 0, 0) and n = (0, 0, 8y, 0) is the separation 
vector of geodesic deviation. Substitution of the curvature expressions (26) 
and (27) into equation (32) yields 

|^N-2BV) 2 N = (33) 

Since the geodesic deviation equation can be considered as a parametrically 
driven harmonic oscillator [9] in the affine parameter A with frequency 



Q(t) = U jRy't'y't' + Ry'x'y'x'(u ) (34) 
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in the turbulent shear layer considered here the frquency is then imaginary 
or Q(X) = iB which leads us to the physical conclusion that the Riemann 
curvature acts locally as a diverging lens in the effective spacetime of the 
GWBH. 



4 Dispersion relation for the turbulent BH 
analog 

By making use of the ansatz 

^(t,f) = V>o(f)e~ <(wt ~*- f) (35) 

and 

0(t, f) = 0o(f)e- i(wot -^ (36) 

and substituting these expressions into the equation (11) one obtains the 
following dispersion relation 



(c 2 k 2 - u 2 ) - [U.V]P -0 2 - 2zu(3 = -^-[q(t,r) - i-^ - 0±] (37) 



where (3 = U.k. From expression (18) we note that there are complex terms 
in frequency ui which implies that we have dissipation terms in the dispersion 
relation. To simplify matters we make use of the ressonance frequency u = ujq 
and ko = k and further assuming that the amplitude of the force <fio is much 
smaller than the amplitute of the potential flow vp one obtains the final 
dispersion relation 



UQ * = fk 2 -[U.V]P - - (38) 



z 2 

Po tpo 

where q is the amplitude of the constante injection of the generator of turbu- 
lence. Note also that there is no dissipation in this dispersion relation. From 
the more general expression (18) one notes that the Lorentz acoustic type 
violation is also present here dueto the turbulent terms. A simple algebraic 
manipulation in these formulas allow us to determine a simple expression for 
the group velocity 

v g = ^ = ±[l-qo]v-iU.k (39) 
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where v := ^- where v p is the phase velocity. The frequency is given by 

w=^[l-ft] (40) 

Here we used the approximation \U\ << v g [10] which was used by Berry et 
al in the investigation of the wavefronts dislocations in the Aharonov-Bohm 
effect in its water wave analog. 

5 conclusions 

A Riemannian shear flow thin layer is shown to yield an effective metric 
for a non-rotating black hole in water turbulent waves. The idea presented 
here extends Schutzhold-Unruh gravity water waves BH. Non-Riemannian 
turbulent flows and their respective BH analogs can also be generated if one 
considers the vorticity of the fluid besides the shear considered here. This 
work is under progress. Another interesting future prospect is the investiga- 
tion of the AB effect around the GWBH as was done around vortices by F. 
Lund [11]. Actually since in these papers dislocation is treated in water a 
Cartan torsion vector could be associated to a non-Riemannian geometry of 
turbulent shear flow like a black hole analog in water waves. In this respect 
knowledge of the non-Riemannian structure of analog models in vortex fluid 
mechanics and acoustics [5] would be very useful. 
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